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University of Helsinki (Finland) 

An accurate determination of the Higgs decay branciiing ratios to bb, cc, gg, WW* 
and r+r~ pairs is important for the study of the Higgs couplings and for determin- 
ing its SM or MSSM nature in the mass range 100 GeV/c^ < Mh < 140 GeV/c^. 
This measurement also represents an important benchmark for the optimisation 
of the detector design. The accuracy on the determination of the Higgs decay 
branching ratio to fermions and WW pairs has been studied using the simulation 
of the detector designed for the TESLA e+e~ linear collider at y's of 350 GeV 
and 500 GeV. The results are discussed in terms of the ability of discriminating 
between the SM and a MSSM neutral Higgs boson and of the predictivity on the 
A^ mass in the MSSM scenario. 



1 Introduction 

The investigation of the production and decay properties of the Higgs boson 
and the determination of its nature are an important part of the physics pro- 
gramme at the e~^e~ Hnear coUider. Furthermore this study sets stringent 
requirements on the response of the experimental apparatus and thus repre- 
sents an ideal benchmark for the optimisation of its design. 

At the linear collider, detection of the SM or MSSM-like neutral Higgs 
boson will be straightforward and the anticipated large statistics and accu- 
rate detector response will enable detailed tests of its production and decay 
characteristicslll. 

This study is Jpased on the simulated response of the detector designed for 
the Tesla CDR0. The track impact parameter resolution has been varied 
to reflect the updates in the Vertex Tracker design El. Signal e+e~ -^ Z^^H^, 
H^i/D and background events have been simulated using the Pythia 5.02 and 
Jetset 7.405 generatorsQ tuned on Lep data for the electroweak and QCD 
variables and on both Lep and Cleo data for heavy flavour decays. Events 
have been generated at y/s of 350 GeV and 500 GeV, corresponding to an 
integrated luminosity oi f L = 100 fb~^. The centre-of-mass energy included 
the smearing effect due to beamstrahlung lI. The results have been scaled to 
the integrated luminosities of 500 fb^^ corresponding to 1 to 2 years (1 year 
= 10^ s) of data taking at the Tesla design luminosity of 3 x 10^"* cm~^ s~^. 
The detector response to charged particle tracks has been studied using a full 



Geant 3.21 simulation and a track fit based on the Kalman filter algorithm. 
The resulting momentum and impact parameter resolutions have been used as 
inputs to a parametric smearing simulation programQ. This program has been 
also used to model the calorimeter response. Events have been reconstructed 
using a dedicated event reconstruction and analysis program. Finally a jet 
flavour tagging algorithm has been applied in order to identify the hadronic 
final states. The reconstruction, analysis and tagging algorithms are based on 
programs developed for the study of the Lep-2 data with the Delphi detector. 

2 Determination of Higgs Branching Ratios 

The accuracy on the Higgs decay branching ratio measursHients at the linear 
collider has been already the subject of several studies IrliJ. With improved 
designs of the Vertex Tracker, more advanced jet fiavour tagging techniques, 
profiting of the experience gained at Lep and Slc, and the large statistics 
available at the Tesla collider these studies move in the domain of precision 
measurements. 

2.1 Event Reconstruction 

This analysis considered neutral Higgs bosons, produced by either the e+e^ — > 
H^Z^ or the e+e" -^ H^uD processes, in the {JJ)h{JJ)z, {JJ)h{^^)z and 
{JJ)hi^v final states, where {J J)h stands both for the two jet hadronic final 
states and for t^t^ and WW* as discussed below. 

After the eventJielection, hadronic jets have been reconstructed using the 
LuCLUS algorithm □. Events have been forced to either two or four jets and 
the jet energies have been rescaled using a 4-C fit by imposing energy and 
momentum conservation at the nominal ^/s. Further selections have been 
applied on the event topology for the three final states under study. Finally the 
compatibility of each reconstructed event with the H'^Z'^, H^vv hypotheses has 
been tested by constructing a x|,j variable based on the reconstructed masses. 
The di-jet Higgs mass resolution has been found to be ajj ~ 2.7 GeV/c^ at 
120 GeV/c^. Candidate Higgs decays have been selected by a cut on this x^^ 
variable. The efficiency of this selection has been estimated to be 25% with 
76% purity, corresponding to a signal to background ratio S/^/B ~ 100, at 
both yft = 350 GeV and 500 GeV. 

The Higgstrahlung cross-section ctzh can be obtained by the study of 
the di-lepton recoil mass, Mu, from Z^ -^ t^i~ . This is necessary for the 
extraction of the absolute Higgs decay branching ratios by flavour tagging 
of the Higgs decay flnal states. In addition, a determination of the Higgs 
production cross-section independent on its decay modes is important in the 



study of possible invisible or exotic Higgs decay channels such as decays to 
SUSY A dedicated study has shown that an accuracy of about 2% on cjzh 
can be obtained for both Mh = 120 GeV/c^ and Mh = 140 GeV/c^ with an 
integrated luminosity of 500 fb^^. 



2.2 iJ° -> hh, cc, 
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Higgs hadronic final states have been selected by imposing hadronic jet quality 
cuts. The separation of the 66, cc^ gg hadronic final states is based on the jet 
flavour tagging response. Due to their relatively long lifetime, cascade 6 — *■ c — *■ 
s decay chain and large invariant mass of beauty hadrons, the decay topology 
with two secondary vertices and the secondary particle invariant mass are very 
distinctive features of b jets. At the linear collider, due to the large boost, 
the 6 hadrons decay on average about 0.8 cm away from the beam collision, 
compared to 0.3 cm at Lep and at the Slc. Furthermore the performances 
of the new generation of Vertex Trackers presently under study for the linear 
collider will allow an efficient reconstruction of secondary charged particle and 
the determination of their point of production. Jet flayour tagging techniques, 
similar to those already successfully applied at LepEJ and at the SlcEj, will 
ensure effective performances due to the favourable kinematics and detector 
response. The tagging algorithm adopted in this analysis combines kinematical 
and topological variables in order to identify bb, cc and light quark di-jets. 
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Figure 1: The distribution of the gg, cc and bb di-jet flavour tag response for Higgs decay 
candidates. The points with error bars show the expected data distribution assuming SM 
couplings for 100 fb^^ and the shaded histograms the response for light quark (medium 
grey), cc (light grey) and bb (dark grey) di-jets normalised to the fitted fractions. 



For each hadronic jet an inclusive secondary vertex search has been per- 
formed using charged particles tracks reconstructed in the Vertex Tracker. 
Neutral particles with large probability of originating from a heavy hadron de- 
cay have been also associated to this secondary vertex. The identification of the 
jet flavour has been based on the following variables: i) probability that all the 
reconstructed charged particles of the jet originated at the event primary ver- 
tex, ii) number of charged particles with an impact parameter larger than twice 
the measurement accuracy, iii) invariant mass and energy of the secondary 
particles, iv) transverse momentum of identified leptons in the jets and v) 
probability that all the secondary particles originated at a common, displaced 
vertex. The likelihood function has been defined as follows: for each of the N 
discriminating variables, the fractions ^^''''(xi), Ff'^{xi) a.ud Ff^ (xi) of 6, candg 
di-jets corresponding to a given value Xi of the i*'' variable, have been extracted 
from samples of Higgs bb, cc and gg decays with equal populations. The bb, 
cc and gg likelihood variables have been computed as the normalised product 
of these individual fractions as: n»=i,jv P^^^ixt)/ Eqq=bb^cc^gg Ut^i^N F^\x^)■ 
Only decays with at least one jet with jet flavour tagging response have been 
considered in this analysis. The fractions of 66, cc and gg Higgs final states 
has been extracted by a binned maximum likelihood fit to these di-jet flavour 
tagging probabilities for the Higgs decay candidates (see Figure |l|) . The back- 
ground has been estimated from simulation over a wide interval around the 
Higgs mass peak and subtracted. It will be possible to study the flavour com- 
position of this background directly on the real data by using the side-bands 
of the Higgs mass peak, thus reducing the possible systematics error from the 
simulation modelling. The fitted fractions agreed well with the input values 
and the performances at y/s ~ 350 GeV and 500 GeV have been found to be 
equivalent. 



Table 1: Accuracy in the determination of the hadronic 
Mh = 120 GeV/c2 and ^ = 350 GeV 



Higgs decay branching ratios for 



^ cTz/f x_Bi?,(ii— i-hadrons) -'/ 

CDR Vtx. I Improved Vtx. 



Channel 



H"/h'' 



bb 



99 



± 0.011 
± 0.134 
± 0.050 



± 0.008 
± 0.080 
± 0.050 



2.3 H° ^ T+T- 

The selection of H^ -^ t^t^ decays started from events not fulfilling the Higgs 
hadronic final state selection. Events have been clusterised in either two or four 
jets with isolated tracks being considered as single particle jets. Similarly to the 

Table 2: Performances of the H ^ tt analysis for Mh = 120 GeV/c^ and y^ = 350 GeV 



Signal H ^ TT 
/lOOfb-i 


Background 
/lOO fb-i 


Eff. 


SBR/BR 
/500 fb-i 


165 


280 


19% 


± 0.057 



case of the jet flavour tagging of hadronic final states, a global tt likelihood has 
been defined by using the response of the following discriminating variables: i) 
number of jets (including isolated particles), ii) number of jets with more than 
2 charged particles in the event, iii) event Thrust, iv) Etot/\/s, v) jet invariant 
mass, vi) number of charged particles per jet and vii) jet impact parameter 
probability. This tt likelihood peaked at one for 7?" -^ t~^t~ decays and at 
zero for background events. H —> tt decays have been selected by a cut on 
this likelihood. The background from other Higgs decay channels has been 
found to be negligeable. Results are summarised in Table 0. 



2.4 H" 



WW* 



The determination of the WW* decay branching ratio has been performed 
using the -ff° -^ WW*^ £vqq channel. This analysis of the H° -^ WW* is 
presented in details in 113. The signal is characterised by the event topology 
and the two-jet recoil mass distribution peaked at the Higgs mass. The results 
are summarised in Table pi 

Table 3: Performances of the H -^ WW* analysis for Mh = 120 GeV/c^ and ,/s = 350 GeV 



H-^WW* 

/lOO fb^i 


Background 

/lOO fb-i 


Eff. 


6BR/BR 
/500 fb-i 


101 


30 


5.3% 


± 0.051 



This measurement, combined with the determination of the Higgstrahlung 
and WW fusicp, cross-sections, can be used to extract the Higgs width with 
good accuracy Ej. 



3 Telling a SM from a MSSM Higgs Boson 

The accuracies on the determination of the decay branching ratio into bb, cc, 
gg, TT and WW*, obtained in this study for the case of a Mh = 120 GeV/c^ 
Higgs boson, are summarised in Table H. Their imphcations on distinguishing 



Table 4: Accuracy in the determination of the Higgs decay branching ratios 



Channel 


5{BR{H -^ X)/BR 




CDR Vtx. 1 Improved Vtx. 


ffV/i" ^ bb 


± 0.024 1 ± 0.024 


HO/h" -^ cc 


± 0.135 1 ± 0.083 


HVh^ - gg 


± 0.055 1 ± 0.055 


H^/h^ _> r+T~ 


± 0.060 


H°/h" -> WW* 


± 0.054 



the SM or MSSM nature of the Higgs boson are discussed in this section. 

The sensitivity of the Higgs decay branching ratios in discriminating be- 
tween tiie SM and MSSM Higgs hypotheses has been already studied in some 
detailsEJ. In MSSM the Higgs decay width to specific final states has the follow- 
ing dependence: T^f ^^^ ex rg:^!is^, T^l^^*'^ ex r?*^22S^ and VM^if^^ 



^wwi with tana — 



- bb cos^ 13 ' 
(M%+M\) sin (i cos (3 
Ml-(M'^ cos2 I3+M\ sin2 /3) 



sin^ 13 



WW 



and tan/3 = ^ being the ratio of 
vacuum expectation values. Therefore deviations in the ratio of branching ra- 



tios such as 



BB.[h~*WW) BB.(h^cc) 



BBjh^gg) 



BRih^bb) ' BRih^bb) ^'^d BRihZbb) f'^O'^ tlieir SM expectations 
could reveal the MSSM nature of the Higgs boson and provide indirect infor- 
mation on the mass of the CP odd A^ Higgs boson. The accuracy of the SM 
predictions has been studied assuming 7nfc(MY(i5)/2) = (4.20 ± 0.11) GeV/c^ 
for the running MS, mb ~ nic ^ (3.40 ± 0.04) GeV/c^, as{mz) = 0.1164 ± 
0.0025 and mtop = (175 ± 0.2) GeV/c^ where the b quaxk mass has been ob- 
tained using a sum rule for the T spectroscopy data cB, the rub^ rric mass 
difference has been derived from Heavy Quark Effective Theory Ej while for 
the top pole mass, mtop, the uncertainty reflects a conservative estimate iQf 
the expected accuracy from the top threshold scan at the linear collider E3. 
Improvements on m(, and mi, — rric, possibly by a factor ~ 2, can be envisaged 
after the study of the data on S-decays from the _B-factories ajid the Lhc. The 
SM Higgs decay branching ratios, including QCD correctionsEJ, have been ob- 
tained using the Hdecay program Ej. Results are summarised in Figure |[ 
The c-quark mass and the as uncertainties limit the predictivity of the cc and 
gg channels to about ±14% and ±7% respectively. On the contrary, the bb, 
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Figure 2: The predicted SM Higgs decay branching ratios shown as the 68% confidence level 
bands with overlayed the measured points using the results of this study. 



T~^T~ and WW* predictions can be obtained with accuracies comparable to, 
or better than, the experimental uncertainties. 

For comparing these SM predictions to those in MSSM, a scan of the 
MSSM parameters phase space has been performed assuming: 

2 < tan/3 < 60 
150 GeV/c2 < Ma < 1100 GeV/c^ 
500 GeV/c2 < MsusY < 1500 GeV/c^ 
-1000 GeV <n< 1000 GeV 
< MlJMg < V6 
0.5 < M~g/MsusY < 1. 
For each set of parameters thej-cesulting /i° mass has been computed using thfi 
diagrammatic two-loop result c3 implemented in the FeynHiggs programed. 
Solutions corresponding to M/jO = (120 ± 2) GeV/c^ have been selected and 
used to compute the /i° decay branching ratios with the Hdecay program, ac- 
counting for squark loops, after forcing M/iO = 120 GeV/c^. The pull quantity 
A{BR) = 1^^ ^ ~f^ — L has been computed for i) BR{h -^ bb)/BR{h -^ 

hadrons), ii) BR{h -> cc)/BR{h -^ hadrons), iii) BR{h -> gg)/BR{h -^ 
hadrons) and iv) BR{h -> bb)/BR{h -^ WW*). From the total x^ the 
SM/MSSM discrimination has been defined as the portion of the Ma — tan /? 
plane with more than 68%, 90% or 95% of the MSSM solutions outside the SM 
95% confidence level region. This study has been repeated for four different 



scenarios: i) 500 fb^^, theory systematics as above with CDR Vertex Tracker, 
ii) 500 fb~^, theory systematics as above with improved Vertex Tracker, iii) 
1000 fb""*^, theory systematics half of the above with CDR Vertex Tracker 
and iv) 1000 fb~^, theory systematics half of the above with improved Vertex 
Tracker. SM/MSSM 90% C.L. discrimination can be achieved up to Ma — 
550 GeV/c^ for the scenario i) improving to Ma — 730 GeV/c^ assuming sce- 
nario iv). Since, with the experimental performances obtained for this analysis, 
the dominant source of uncertainties is due to the theory systematics, there 
is not a significant difference in the excluded region between the two assumed 
performances for the Vertex Tracker. The results are exemplified in Figure ^ 
From this comparison it follows that the sensitivity to the SM/MSSM nature of 
a neutral Higgs boson is mainly limited by the theory uncertainties rather than 
by the envisaged detector response and the experimental accuracies obtained 
in the present study. 

TESLA L = 500 fb"" 
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Figure 3: The SM/MSSM discrimination region in the Ma — tan/3 plane for Mh = 
120 GeV/c^ and 500 fb"!. The lines define the regions with 68% (light grey), 90% (medium 
grey) and 95% (dark grey) of the MSSM solutions distinguishable at the 95% confidence 
level from the SM H^ boson. 



It is also possible to use the accurate determination of the Higgs decay 
branching ratios for an indirect estimate of the mass Af^o in the MSSM case. 
Sensitivity to the A^ mass arises from the MSSM corrections to the Higgs 
couplings discussed above and it is of special interests for those masses above 
the kinematical limit for direct e+e^ -^ hPA^^ H^A^ production E3. The 
analysis has been performed assuming given sets of measured values for the 



BR(/i -^ cc + gg)/BR{h -> bb) and BR{h -^ WW*)/BR{h -^ bb) ratios. The 
A^ mass has been varied together with the other MSSM parameters within 
the range compatible with the measured branching ratios allowing for their 
total uncertainty. The range of values of Ma for the accepted MSSM solu- 
tions corresponded to an accuracy of 70 GeV/c^ to 90 GeV/c^ for the indirect 
determination of Ma in the mass range 300 GeV/c^ < Ma < 500 GeV/c^. 

4 Conclusions 

The accuracy on the determination of the decay branching ratios for a neutral 
Higgs boson with Mh« = 120 GeV/c^ has been studied at a y/s = 350 GeV 
and 500 GeV e+e^ linear collider. The measured branching ratios allow to 
determine the Higgs decay width with good accuracy, to distinguish a MSSM 
h'^ boson from the SM iJ" boson and to indirectly determine the A" mass up 
to 600-700 GeV/c^. Significant improvements in the knowledge of the b and c 
quark masses and good control of QCD corrections are necessary in order to 
preserve the predictive potential of the cc and gg decays for Mh < 120 GcV/c^. 
For higher Higgs masses, the combination of the bb and WW* decay branching 
ratios provides the bulk of the information on the neutral Higgs nature and 
the A'^ mass. 
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